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Abstract 

Generalizing ideas of MacKay, and MacKay and Saffman, a necessary condition for the 
presence of high-frequency (he., not modulational) instabilities of small-amplitude periodic 
solutions of Hamiltonian partial differential equations is presented, entirely in terms of the 
Hamiltonian of the linearized problem. With the exception of a Krein signature calculation, 
the theory is completely phrased in terms of the dispersion relation of the linear problem. 
The general theory changes as the Poisson structure of the Hamiltonian partial differential 
equation is changed. Two important cases of such Poisson structures are worked out in full 
generality. An example not fitting these two important cases is presented as well, using a 
candidate Boussinesq-Whitham equation. 


1 Introduction 


It is expected that much of the dynamics of the small-amplitude solutions of a partial differential 
equation (PDE), including their stability or instability, is dictated by the study of a linearized 
(about a trivial solution, say rt = 0) problem. In this article, we focus specifically on the spectral 
stability of periodic traveling-wave solutions of Hamiltonian PDEs as they bifurcate away from 
a trivial solution. Our work follows earlier ideas of MacKay [35] and MacKay and Saffman [36] . 
We start from an autonomous Hamiltonian system of PDEs [2], f.e.. 


SH 

ut = J^- 
ou 


( 1 ) 


Here and throughout, indices involving x or t denote partial derivatives. Eurther, u = 
{ui{x,t),... is an M-dimensional vector function defined in a suitable function 

space, and J is a Poisson operator mm- More details and examples are given below. Ei- 
nally, H = Pf (tt, Ux, ■ ■ ■)dx is the Hamiltonian, whose density % depends on u and its spatial 
derivatives, defined for x G D. We consider only the stability of periodic solutions, thus D is 
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Figure 1: A cartoon of the bifurcation structure of the traveling waves for a third-order (M = 3) 
system: solution branches bifurcate away from the trivial zero-amplitude solution at specihc 
values of the traveling wave speed c. 


any interval of length L, the period. Note that for some of our examples % will depend on 
spatial derivatives of u of arbitrary order. 

To investigate the stability of traveling wave solutions of this system, we reformulate 0 
in a frame moving with speed c, using the transformation x = x — ct, t = t, and considering 
solutions u{x,t) = U{x) (successively omitting hats). This leads to 


_ 

CUx — r r 5 

OU OU 


( 2 ) 


for a modified Hamiltonian He- Traveling wave solutions are solutions of the ordinary differential 
system 


r. JH ^ JHc 
-eUx-J^ ^ 


(3) 


Thus if J is invertible, traveling waves are stationary points of the Hamiltonian He- The sys¬ 
tem Q typically has the zero (trivial) solution for a range of c values. The small-amplitude 
solutions whose stability we investigate bifurcate away from these trivial zero-amplitude solu¬ 
tions at special values of the speed parameter c, as is schematically shown in Fig. [T] It is our 
goal to see to what extent anything can be said about the stability of the small-amplitude so¬ 
lutions (with amplitudes in the shaded regions of Fig.[^ from knowledge of the zero-amplitude 
solutions at the bifurcation point. An outline of the steps in this process is as follows. 


1. Quadratic Hamiltonian. A linear system of equations is obtained by linearizing the 
system 0 around the zero solution: let u = ev + o(e) and omit terms of order o(e). 
Alternatively, if J is independent of u and its spatial derivatives, one may expand the 
Hamiltonian He as a function of e and retain its quadratic terms. The resulting Hamilto¬ 
nian He of the linearized system is the starting point for the next steps. 

2. Dispersion relation. The linearized system has constant coefficients and is easily solved 
using Fourier analysis. The dispersion relation F{uj, k) = 0 governs the time dependence 
of the solutions. It is obtained by investigating solutions whose spatial and temporal 
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dependence is proportional to exp(iA:x — itot). Here F{u},k) = 0 is of degree M in cj. It 
is a fundamental assumption of our approach that all solutions ujj{k) {j = 1, ... ,M) of 
F(uj, k) = 0 are real for A: G M. The dispersion relation can be expressed entirely in terms 
of the coefficients appearing in the quadratic Hamiltonian For periodic systems of 
period L, the values of k are restricted to be of the form 2ttN/L, iV G Z. 

3. Bifurcation branches. The values of the phase speed cj = ojj/k for which nontrivial 
solutions bifurcate away from the zero-amplitude solution are determined by the condition 
that the zero solution is not the unique solution to the Fourier transformed problem. In 
effect, this is the classical bifurcation condition that a Jacobian is singular. This simple 
calculation determines the bifurcation branch starting points explicitly in terms of the 
different solutions to the dispersion relation. In what follows, we follow the first branch, 
starting at ci, without loss of generality. 

It is assumed that only a single non-trivial bifurcation branch emanates from a bifurcation 
point. Although more general cases can be incorporated, we do not consider them here. 
Further, we fix the period of the solutions on the bifurcation branch (usually to 27r). Other 
choices can be made. Instead of varying the amplitude as a function of the speed for fixed 
period, one could hx the speed and vary the period, etc. The methods presented can be 
redone for those scenarios in a straightforward fashion. 

4. Stability spectrum. The spectrum of the linear operator determining the spectral 
stability of the zero solution at the bifurcation point on the first branch is calculated. 
Since this spectral problem has constant coefficients, this calculation can be done explicitly. 
Again, this is done entirely in terms of the dispersion relation of the problem. Using a 
Floquet decomposition (see [El El]), the spectrum is obtained as a collection of point 
spectra, parameterized by the Floquet exponent fj. G {—'KjL,'KjL\. Due to the reality of 
the branches of the dispersion relation, the spectrum is confined to the imaginary axis. 
In other words, the zero-amplitude solutions are spectrally stable. The use of the Floquet 
decomposition allows for the inclusion of perturbations that are not necessarily periodic 
with period L. Instead, the perturbations may be quasiperiodic with two incommensurate 
periods, subharmonic (periodic, but with period an integer multiple of L), or spatially 
localized [El EH El]. 

5. Collision condition. Given the explicit expression for individual eigenvalues A, it is easy 
to find the conditions for which eigenvalues corresponding to different parameters (Floquet 
exponent, branch number of the dispersion relation, etc.) coincide on the imaginary axis. 
This is referred to as the collision condition. Once again, it is given entirely in terms of 
the dispersion relation. 

It is a consequence of the Floquet theorem [E] that collisions need to be considered only 
for spectral elements corresponding to the same value of the Floquet exponent since the 
subspaces of eigenfunctions for a fixed Floquet exponent are invariant under the flow of 
the linearized equation. 

6. Krein signature. Having obtained the stability spectrum at the starting point of the 
bifurcation branches, we wish to know how the spectrum evolves as we move up a bifur¬ 
cation branch. One tool to investigate this is the Krein signature |32l |33l IMl ESI EE]- In 
essence, the Krein signature of an eigenvalue is the sign of the Hamiltonian of the lin¬ 
earized system evaluated on the eigenspace of the eigenvalue. Different characterizations 
are given below. If two imaginary eigenvalues of the same signature collide as a parame¬ 
ter changes, their collision does not result in them leaving the imaginary axis. Thus the 
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Figure 2: Colliding eigenvalues in the complex plane as a parameter is increased. On the left, 
two eigenvalues are moving towards each other on the positive imaginary axis, accompanied 
by a complex conjugate pair on the negative imaginary axis. In the middle, the eigenvalues 
in each pair have collided. On the right, a Hamiltonian Hopf bifurcation occurs: the collided 
eigenvalues separate, leaving the imaginary axis (implying that the two Krein signatures were 
different). 


collision of such eigenvalues does not result in the creation of unstable modes. In other 
words, it is a necessary condition for collisions to lead to instability that the Krein signa¬ 
ture of the colliding eigenvalues is different. This scenario is illustrated in Fig. That 
figure also illustrates the quadrifold symmetry of the stability spectrum of the solution of 
a Hamiltonian system: for each eigenvalue A G C, A*, —A and —A* are also eigenvalues. 
Here A* denotes the complex conjugate of A. It should be noted that the occurrence of 
a collision is required for eigenvalues to leave the imaginary axis, due to the quadrifold 
symmetry of the spectrum. 

Thus we calculate the Krein signature of any coinciding eigenvalues, obtained in Step 5. If 
these Krein signatures are equal, the eigenvalues will remain on the imaginary axis as the 
amplitude is increased. Otherwise, the eigenvalues may leave the imaginary axis, through 
a so-called Hamiltonian Hopf bifurcation |45j , resulting in instability. Thus we establish a 
necessary condition for the instability of periodic solutions of small amplitude. The Krein 
signature condition cannot be expressed entirely in terms of the dispersion relation, and 
the coefficients of are required as well. Please refer to the next two sections for details. 

Although all calculations are done for the zero-amplitude solutions at the starting point of 
a bifurcation branch, the continuous dependence of the stability spectrum on the parameters in 
the problem [2S], including the velocity of the traveling wave or the amplitude of the solutions, 
guarantees that the stability conclusions obtained persist for solutions of small amplitude. Thus 
meaningful conclusions about solutions in the shaded regions of the bifurcation branches of Fig.[^ 
are reached, despite the fact that we are unable to say anything about the size of the maximal 
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amplitude for which these conclusions are valid. 

Remarks. 

• Throughout this introduction, and in fact throughout much of this manuscript, our em¬ 
phasis is on generality and usability as opposed to rigor. As a consequence some of the 
statements made above are necessarily vague: more precise statements would limit the 
generality aimed for. Within the context of more specific examples, more precision may be 
possible. Along the same lines, we have limited ourselves to the most generic case of two 
eigenvalues colliding on the imaginary axis. Lastly, since we are assuming an initial (ie., 
for a starting value of 0 for the amplitude parameter) situation that is neutrally stable, 
our only interest is in collisions on the imaginary axis. If eigenvalues are present off the 
imaginary axis for the zero-amplitude solution, such solutions are already spectrally un¬ 
stable, and the continuous dependence of the spectrum on its parameters [26] guarantees 
instability for solutions of small amplitude as a consequence. 

• If eigenvalues collide at the origin, their Krein signature is zero (this follows immedi¬ 
ately from its definition) and no conclusion can be drawn about whether or not colliding 
eigenvalues leave the imaginary axis or not. Because of this, the methods discussed here 
allow only for the study of so-called high-frequency instabilities. Indeed, if eigenvalues 
collide away from the origin, the unstable perturbations do not only grow in magnitude 
exponentially, they also display an oscillation of frequency u) that is equal to the non-zero 
imaginary part of the colliding eigenvalues. In this sense high frequency is more accu¬ 
rately described as non-zero frequency. We use the high-frequency name to distinguish 
the instabilities investigated here from the modulational instability |5l|, which is a con¬ 
sequence of eigenvalues colliding at the origin. Further, as seen in the examples below, 
often the presence of one high-frequency instability is accompanied by a sequence of such 
instabilities with increasing frequencies tending to infinity. 

The study of the modulational instability requires a different set of techniques, see for 
instance miiziisii for different classical approaches. More recently, a general framework 
was developed by Bronski, Johnson and others (see e.g. |9l|27|). Collisions at the origin 
are common, since any Lie symmetry of the underlying problem gives rise to a zero 
eigenvalue in the stability problem |23l[3l|. For this same reason, such collisions typically 
involve more than two eigenvalues (for instance, the original Benjamin-Feir instability 
involves four mn), which is one of the reasons why their treatment is often complicated: 
its analysis is not only technical but also involves tedious calculations, as exemplified 
recently in a paper by Hur and Johnson m on the modulational instability for periodic 
solutions of the Whitham equation. 

• Following the Floquet decomposition, we use eigenfunctions given by a single Fourier 
mode since the linear stability problem for the zero-amplitude solutions has constant co¬ 
efficients. Lastly, the calculation of the Krein signature involves only the finite-dimensional 
eigenspace of the eigenvalue under consideration. Thus, in essence, all the calculations 
done in this paper are finite dimensional, as they are in [36|, for instance. As was noted 
there, for equations with real-valued solutions, it is necessary to consider the eigenfunc¬ 
tions corresponding to Floquet exponent /r and —/r simultaneously, since the eigenfunctions 
with —jj are the complex conjugates of those with /x. 

• We call an equation dispersive if all branches of its dispersion relation (jj[k) are real valued 
for A: G M. It is easy to see that there exist linear, constant-coefficient Hamiltonian systems 
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that are not dispersive. An explicit example is 

] Q.t — Qxxi /4\ 

\ Pt = -Pxx- ^ 

This (admittedly bizarre) system is Hamiltonian with canonical Poisson structure 

^=(-1 o )’ ® 

and Hamiltonian H = — q^Pxdx. The branches of its dispersion relation are given by 

wi,2(A:) = Tik'^- 

It is an interesting question whether there exist dispersive systems that are not Hamilto¬ 
nian. One may be tempted to consider an example like 


( qt = aqx, 
\ Pt = bpx. 


( 6 ) 


For a, 6 G M this system is dispersive (cui = —ak, UJ 2 = —bk), but it is not Hamiltonian 
with canonical structure. However, allowing for the noncanonical structure 


J = 


dx 0 \ 
0 d, 


(7) 


the system is easily found to have the Hamiltonian H = ^ {aq^+bp^)dx. We conjecture 
that no systems exist that are dispersive but not Hamiltonian. As demonstrated by the 
example above, the question is difficult to analyze, since different forms of the Poisson 
operator have to be allowed. 


• The investigation of the Hamiltonian structure of the zero-amplitude solution follows ear¬ 
lier work by Zakharov, Kuznetsov, and others. A review is available in [53]. The stability 
of the trivial solution is also investigated there, with stability criteria given entirely in 
terms of the branches of the dispersion relation. Using canonical perturbation theory, 
Hamiltonians containing higher-than-quadratic terms are considered. This is used to con¬ 
sider the impact of nonlinear effects on the dynamics of the trivial solution. At this point, 
connections to resonant interaction theory 131231101 become apparent, as they will in 
what follows. Although the physical reasoning leading to resonant interaction theory and 
that leading to the criteria established below is different, it is clear that connections exist. 


2 Motivating example 

Our investigations began with the study of the so-called Whitham equation |39| , [53 page 368]. 
The equation is usually posed on the whole line, for which the equation satisfied by u{x, t) is 

/ OO 

K{x - y)uy{y,t)dy = 0. (8) 

-OO 

The term N(u) denotes the collection of all nonlinear terms in the equation. It is assumed that 
lim^-^o = 0. The last term encodes the dispersion relation of the linearized Whitham 

equation. The kernel K{x) is the inverse Fourier transform of the phase speed c{k), where 
c{k) = Lo{k)/k, with uj{k) the dispersion relation. Here c{k) is assumed to be real valued and 
nonsingular for /c G M. Thus 
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1 

K{x) = i c{k)e^^^dk, (9) 

27r7-oo 

where f denotes the principal value integral. Depending on c(k), this equation may have to be 
interpreted in a distributions sense [20lll3]. Letting u ~ exp{ikx — uj{k)t) , it is a straightforward 
calculation to see that the dispersion relation of the linear Whitham equation is u!{k). In fact, 
the linear Whitham equation is easily seen to be a rewrite of the linear evolution equation [1] 

ut = -iu}{-idx)u, ( 10 ) 

where oj(—idx) is a linear operator with odd symbol uj{k): oj{k) = —uj{—k), the dispersion 
relation considered. Indeed, letting uj{—idx) act on 

1 

u{x,t) = — / e^^^u{k,t)dk, ( 11 ) 

2vr 7_oo 

and replacing u{k, t) by 

/ OO 

e~"''^u{y,t)dy, ( 12 ) 

-OO 

the linear part of Q is obtained after one integration by parts. We restrict our considerations 
to odd dispersion relations, to ensure the reality of the Whitham equation. 

One of Whitham’s reasons for writing down the Whitham equation [491150] was to describe 
waves in shallow water (leading to the inclusion of a KdV-type nonlinearity N(u) ~ uux) that 
feel the full dispersive response of the one-dimensional water wave problem (without surface 
tension), for which 


uj‘^{k) = gktanh{kh), 


(13) 


with g the acceleration of gravity and water depth h. It is common to choose c{k) = uji(k)/k > 0 


m 


so that 0 Ji{k) is the root of (13) with the same sign as k. In what follows, we refer to 


this choice as the Whitham equation. The stability of periodic traveling wave solutions of the 
Whitham equation has received some attention recently. Notable are [19] . where the focus 
is on solitary waves, and where the modulational instability of small-amplitude periodic 
solutions is emphasized, although the main result for the high-frequency instabilities discussed 
in Example 3.2 is included as well. Most recently, the spectral stability of periodic solutions 


of the Whitham equation was examined in mi- The goal of considering an equation like ^ 
as opposed to the Korteweg-de Vries equation or other simpler models is to capture as much 
of the dynamics of the full water wave problem as possible, without having to cope with the 
main difficulties imparted by the Euler water wave problem [46] {e.g., it is a nonlinear free 
boundary-value problem, the computation of its traveling-wave solutions is a nontrivial task, 
etc). One of the important aspects of the dynamics of a nonlinear problem is the (in)stability of 
its traveling wave solutions. It was shown explicitly in m that periodic traveling wave solutions 
of the one-dimensional Euler water wave problem are spectrally unstable for all possible values 
of their parameters h, g, amplitude, and wave period. The nature of the instabilities depends 
on the value of these parameters. As is well known, waves in deep water are susceptible to the 
Benjamin-Eeir or modulational instability (see |5Tj for a review). In addition, waves in both 
deep and shallow water of all non-zero amplitudes are unstable with respect to high-frequency 
perturbations: these are perturbations whose growth rates do not have a small imaginary part, 
resulting in oscillatory behavior in time, independent of the spatial behavior of the unstable 
modes. The work in m does not reveal any high-frequency instabilities for solutions of small 
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amplitude in water that is shallow, in the context of the Whitham equation. Thus an important 
aspect of the Euler water wave dynamics is absent from Q. We will provide an analytical 
indication that the Whitham equation misses the presence of these instabilities, while explaining 
why they are missed. This explanation leads to a way to address this problem. 

For suitable N{u), the Whitham equation Q is a Hamiltonian system. In fact, for our 
considerations, it matters only that the linearized Whitham equation is Hamiltonian. The 


Lagrangian structure with the dispersion relation given by wi(fe) as in (13) was already written 
down by Whitham in |48j . from which the Hamiltonian structure easily follows. Explicitly, for 
the linearized Whitham equation posed on the whole line, for any odd uj{k) we have 


H = - 


/‘OO poo 

- / K{x - y)u{x)u{y)dxdy, 

^ J —OO j —OO 


(14) 


with J = dr- Then 


ut = d, 


5H_ 

Su 


(15) 


If instead the linearized equation is posed with periodic boundary conditions u{x+L, t) = u{x, t), 

(16) 


it follows immediately from (10) that we have 

fL/2 


Ut 


+ K{x-y)u{y)dy = 0, 
J-L/2 


where we have used a Fourier series instead of a Fourier transform. Further, 


^ OO 

K{x) = - ^ c(%)e* 


(17) 


j=-oo 


and kj = 27rj/L, j G Z. The Hamiltonian formulation for the periodic Whitham equation (16) 
is also given by (15), but with 


H = -l 


fL/2 rL/2 

/ / K{x — y)u{x)u{y)dxdy. (18) 

'-L/2 J-L/2 

In fact, a formal L —)• oo immediately results in the recovery of the equation posed on the whole 
line. Thus the Whitham equation and its periodic solutions fit in to the framework developed 
in this manuscript. It is one of many examples we use below. Other notable examples are the 
Euler water wave problem (as expected, allowing us to check our results with those of MacKay &: 
Saffman [3^), the KdV equation, the Sine-Gordon equation, etc. We are particularly interested 
in the comparison between the results for the Euler water wave problem and those for the 
Whitham equation. 


The results of Examples 3.2 and 4.2 show that the Whitham equation cannot possess the 
high-frequency instabilities present in the water wave problem. This leads us to propose a new 
model equation, a so-called Boussinesq-Whitham or bidirectional Whitham equation. This equa¬ 
tion is shown to at least satisfy the same necessary condition for the presence of high-frequency 
instabilities as the water wave problem, and these high-frequency instabilities originate from 
the same points on the imaginary axes as they do for the Euler equations. However, it is easily 
seen that this equation is ill posed for solutions that do not have zero average. As such it is 
a poor candidate to replace the Whitham equation ([^ as a shallow-water equation with the 
correct dispersive behavior. 


Remark. If oj{k) is not odd, then the function u solving the linear Whitham equation (10) 
is necessarily complex. The problem is Hamiltonian: 
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(19) 


ut = -i 


. 5H 
6u*' 


where * denotes the complex conjugate. The Poisson structure is 


J = 


0 -i 
i 0 


( 20 ) 


although usually one writes only the first of the two evolution equations, omitting the equation 
that is the complex conjugate of (19). The Hamiltonian is given by 


f27r 


H = 


u*uj{—idx)udx. 


( 21 ) 


A real formulation in terms of the real and imaginary parts of u is possible as well (using 
the canonical transformation u = {q + ip)/^/2), resulting in a canonical Hamiltonian struc¬ 
ture. The linear Whitham equation (19) is often rewritten in Fourier space using the discrete 
Fourier transform, due to the periodic boundary conditions. This leads to the Hamiltonian 


3 Scalar Hamiltonian PDEs 

In this section, we investigate the stability of 27r-periodic traveling wave solutions of Hamiltonian 
systems of the form 


ut = dx^, ( 22 ) 

du 

where u{x,t) is a scalar real-valued function. Thus J = dx- Since this Poisson operator is 
singular, all equations of this form conserve the quantity udx, which is the Casimir for this 
Poisson operator. Systems of this form include the Korteweg- de Vries equation [22l|52] and its 
many generalizations, the Whitham equation ([^, and many others. As mentioned above, our 
only interest is in the linearization of these equations around their trivial solution. We write 
the quadratic part of H as 


^ 


c27r 

, V 

J n=0 


CkfiUx^xdx , 


(23) 


where the coefficients an G M. As before, indices on u denote partial derivatives. Specifically, 
Unx denotes d'^u. 

For most examples, the number of terms in (23) is finite, and all but a few of the coeffi¬ 


cients an are nonzero. For the Whitham equation J8), the number of nonzero terms is infinite, 
but convergence is easily established. Note that (1^) is the most general form of a quadratic 
Hamiltonian depending on a single function. Indeed, a term in containing UmxUnx, with m 
and n positive integers, may be reduced to a squared term using integration by parts. 


Using the notation (23), the linearized equation is 


^ ^ ( 1 ) CXn'^^(2n+l)x- 

n=0 


(24) 


We proceed with the six steps outlined in the introduction. 
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1. Quadratic Hamiltonian. The modified Hamiltonian is given by 



( 25 ) 


2 . 


Dispersion relation. For 

single branch: 


equations of the form (22), the dispersion relation has only a 




(26) 


n=0 


The absence of even powers of A: in (26) is due to our imposition that (24) is a conservative 


equation, i.e., there is no dissipation. All integers are allowable k values, since we have 


equated the period to be 27r. The equation (24) may be written as 


ut = —iuj{—idx)u. 


(27) 


3. Bifurcation branches. Since (22) is scalar, only one branch can bifurcate away from 


the trivial solution. To find the corresponding value of c, we write (24) in a moving frame 


as 


Ut — cux = ioj{idx)u. (28) 

This equation has its own dispersion relation given by 

Vt{k) = u:{k) — ck, (29) 


obtained by looking for solutions of the form u = ex.p{ikx — fOt). Letting u = 
Yl’^=-oo exp(inx)un, it follows that dtUn = —i^{n)un- Thus a nonzero stationary so¬ 
lution may exist provided Q{N) = 0, for N € N, N ^ 0. We have used the oddness of 
Q{N) to restrict to strictly positive values of N. Thus the starting point of the bifurcation 
branch in the (speed, amplitude)-plane is (c, 0), where c is determined by 


c = 


u{N) 

N 


(30) 


for any integer > 0. Choosing iV > 1 implies that the fundamental period of the solu¬ 
tions is not 27r, but 27r/A^. In practice, we choose A" = 1. A Fourier series approximation 
to the explicit form of the small-amplitude solutions corresponding to this bifurcation 
branch may be obtained using a standard Stokes expansion [ill [50] . 


4. Stability spectrum. In order to compute the stability spectrum associated with 
the zero-amplitude solution at the start of the bifurcation branch, we let u{x, t) = 
I/(x) exp(At)-|-c.c., where c.c. denotes the complex conjugate of the preceding term. As 
usual, if any A are found for which the real part is positive, the solution is spectrally 
unstable m- All bounded eigenfunctions U (x) may be represented as 


OO 

U{x)= Y. 

n=—oo 


( 31 ) 
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where fi G (—1/2,1/2] is the Floquet exponent. Such a representation for U{x) is valid 
even for solutions on the bifurcation branch of nonzero amplitude m- Since (|2^ is a 


problem with constant coefficients, only a single term in (31) is required. We obtain 


= —iVt{n + ix) = —iuj{n + n) + i{n + /r)c, n G Z. 


(32) 


As expected, all eigenvalues are imaginary and the zero-amplitude solution is neutrally 
stable. For a fixed value of /r, ( |32[ ) gives a point set on the imaginary axis in the complex 
A plane. As ^ is varied in (—1/2,1/2], these points trace out intervals on the imaginary 
axis. Depending on uj{k), these intervals may cover the imaginary axis. 


5. Collision condition. The most generic scenarios for two eigenvalues given by (32) to 
collide are that (i) two of them are zero, and they collide at the origin, and (ii) two of 
them are equal, but nonzero. We ignore the first possibility, since the next step proves 
to be inconclusive for this case, as discussed in the introduction. The second possibility 


requires = Xm\ for some m,n G Z, m / n, fixed ^ G (—1/2,1/2], and xi^^ 
This may be rewritten as 


Xt^ 0 . 


uj(n + fx) — uj(m +/x) co{N) 


n — m 


N 


m, n G Z, m / n and /i G (—1/2,1/2]. (33) 


This equation has an elegant graphical interpretation: the right-hand side is fixed by the 
choice of N, fixing the bifurcation branch in Step 3. It represents the slope of a line 
through the origin and the point {N,uj{N)) in the {k,uj) plane. The left hand side is 
the slope of a line in the same plane passing through the points (n -|- fx,Lj{n + fx)) and 
(m -|- /X, uj{m + fx)), see Fig. 

Even though the graph of the dispersion relation admits parallel secant lines, this is 
not sufficient for a solution of (33), as it is required that their abscissas are an integer 


apart. Nevertheless, the graphical interpretation can provide good intuition for solving 
the collision condition, which typically has to be done numerically. 


6. Krein signature. The Krein signature of an eigenvalue is the sign of the Hamiltonian 
evaluated on the eigenspace associated with the eigenvalue. We are considering two simple 
eigenvalues colliding, thus the eigenspace for each eigenvalue consists of multiples of the 
eigenfunction only. To allow for eigenfunctions of the form an exp(i(n -|- fx)x + xli^H)+c.c, 


which are not 27r-periodic (unless /x = 0), it is necessary to replace the integral in (23) 


with a whole-line average. More details on this process are found, for instance, in HZ]. A 
simple calculation shows that the contribution to from the (n, /x) mode is proportional 
to the ratio of H(n -I- /x)/(n -|- /x): 




2 H(n -|- /x) 
n + IX 


(34) 


Other terms are present in the Hamiltonian density, but they have zero average. The sign 
of this expression is the Krein signature of the eigenvalue An^^. Thus a necessary condition 
for aI^^ and Xm^ to leave the imaginary axis for solutions of non-zero amplitude is that 


the signs of (34) with {n,ix) and (m,^) are different, contingent on /x, m and n satisfying 


(33). Explicitly, this condition is 


sign 


oo{n + ix) 
n + fx 


— c 


7 ^ sign 


co{m + fx) 
m + fx 


— c 


(35) 
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Figure 3: The graphical interpretation of the collision condition (33). The solid curve is the 
graph of the dispersion relation uj{k). The slope of the dashed line in the first quadrant is the 
right-hand side in (33). The slope of the parallel dotted line is its left-hand side. 


Alternatively, the product of the left-hand side and the right-hand side should be negative. 
Using (33), (35) becomes 


(n -I- n){m -I- ^) < 0, 

or, provided mn ^ 0, and using that /x G (—1/2,1/2], 

nm < 0. 


(36) 


(37) 


Remarks. 


It is clear from (34) why our methods do not lead to any conclusions about collisions of 
eigenvalues at the origin. If = 0, then fl(n + fi) = 0, and the contribution to the 
Hamiltonian of such a mode vanishes. As a consequence, the associated Krein signature 
is zero. 


• When the theory of |35] is restricted to the case of solutions of zero-amplitude, so as to 
recover the constant coefficient stability problem, the graphical stability criterion given 
there coincides with the one presented here. 


We conclude our general considerations of this section with the following summary. 


Assume that the linearization of the scalar Hamiltonian system (22) is dispersive (i.e., its 
dispersion relation u!(k) is real valued for k G M). Let be a strictly positive integer. Consider 
27r/iV-periodic traveling wave solutions of this system of sufficiently small-amplitude and with 
velocity sufficiently close to u}{N)/N. In order for these solutions to be spectrally unstable 
with respect to high-frequency instabilities as a consequence of two-eigenvalue collisions, it is 
necessary that there exist n, m € Z, n ^ m, fj, € (—1/2,1/2] for which 
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( 38 ) 


such that 


and 


uj{n + p) uj{N) 
n+n ^ N ’ 


u{m + /i) uj{N) 
m + p, N 


uj{n + ^) — uj{m + ^) uj{N) 
n — m N ' 


(39) 


(m +/u)(n +/x) < 0. (40) 

Next we proceed with some examples. 

3.1 The (generalized) Korteweg - de Vries equation 

We consider the generalized KdV (gKdV) equation 

Uf (TU Ux Uxxx — 0 ; (^ 1 ) 

where we restrict n to integers 1 or greater. Here cr is a constant coefficient, chosen as conve¬ 
nient. Important special cases discussed below are the KdV equation (n = 1) and the modified 
KdV (mKdV) equation (n = 2). Many of the details below extend easily to more general 
nonlinearities, with the main requirement being that the linearized equation is ut + Uxxx = 0. 
The stability of periodic solutions of the gKdV equation has received some attention recently 
inmniEii^E]. For the integrable cases n = 1 and n = 2, more detailed analysis is possible, see 
laiiaiie]. We do not claim to add anything new to these discussions, but we wish to use this 
example to illustrate how the six-step process outlined in this section leads to easy conclusions 
before moving on to more complicated settings. 


1. The modified Hamiltonian is given by 


1 

He = 2 Jq cu^)dx. (42) 

2. The dispersion relation is 


w = -k\ 

3. Bifurcation branches in the (c, amplitude)-plane start at (c, 0), with 


(43) 


io{k) 

~k~ 


= -k\ 


(44) 


Since we desire 27r periodic solutions, we choose k = 1. Any choice k = N, where V is a 
non-zero integer is allowed. Choosing k = 1, bifurcation branches start at (—1,0). 

For the integrable cases n = 1 (KdV) and n = 2 (mKdV), these bifurcation branches may 
be calculated in closed form. For the KdV equation in a frame traveling with speed c, the 
27r-periodic solutions are given by (with cr = 1) 


12k?K^(k) 2 ( K{k)x 

u =-^-^cm - 

TT^ V 



(45) 
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Figure 4: The amplitude vs. c bifurcation plots for the traveling-wave solutions of the gener¬ 
alized KdV equation (41). (a) The KdV equation, n = 1, for the cnoidal wave solutions (45). 
(b) The mKdV equation, n = 2, for the cnoidal wave solutions (47). Lastly, (c) shows the 
bifurcation plot for the snoidal wave solutions (48) of mKdV, n = 2. Note that all bifurcation 
branches start at (—1,0), as stated above. 


where cn denotes the Jacobian elliptic cosine function and K{k) is the complete elliptic 
integral of the first kind [T8l|39]. Further, 

c(^) = i^(2K2-l), (46) 

TT^ 

resulting in an explicit bifurcation curve {c{k) , 12 k‘^K‘^(k) / 'k‘^) , parameterized by the el¬ 
liptic modulus K G [0,1). This bifurcation curve is shown in Fig. [^. 

For the mKdV equation (n = 2), different families of traveling-wave solutions exist |16j . 
We consider two of the simplest. For <7 = 3 (focusing mKdV), a family of 27r-periodic 
solutions is given by 


2\/2kK{k) /2K{k)x \ 

u = -cn -, K , 

TT \ TT J 


(47) 


with c(k) given by (46), resulting in an explicit bifurcation curve (c(k), 2\/2KiL(K)/7r), 
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(a) 


(b) 


Figure 5: (a) The imaginary part of G (—0.7,0.7) as a function of /r G [—1/4,1/4). (b) 

The curves Q{k + n), for various (integer) values of n, illustrating that collisions occur at the 
origin only. 


parameterized by the elliptic modulus k G [0,1). This bifurcation curve is shown in 
Fig. 1)3. It should be noted that a solution branch exists where the solution is expressed 
in terms of the Jacobian dn function [T8l[39]: u = ^/{2)K{K)dn{K{^^)x/^^,K)/7^, but this 
solution does not have a small-amplitude limit and our methods do not apply directly to 
it. Rather, the solutions limit to the constant solution u = l/\/2 as k —)• 0. A simple 
transformation v = u — 1 /\/2 transforms the problem to one where our methods apply. 

For a = —3 (defocusing mKdV), a period 27r solution family is 


2\/2kK(k) f 2 K(k)x a 
u = -sn -, K , 

TT V ^ / 


(48) 


with c{k) = —4(1 -|- Here sn is the Jacobian elliptic sine function [181139| . 

resulting in an explicit bifurcation curve (c(k), 2y/2 KK{K)/tt), parameterized by the elliptic 
modulus K G [0,1). This bifurcation curve is shown in Fig. |^. 


4. The stability spectrum is given by (32), with uj{k) = —k^ and c = —1, resulting in 


= i{n + n){l + {n + m)^)- 


(49) 


These eigenvalues cover the imaginary axis, as n and /i are varied. The imaginary part of 
this expression is displayed in Fig. [^. For the sake of comparison with Fig. 2 in [6] we let 
[I G [—1/4,1/4), which implies that n is any half integer. The results of Fig. 2 in [6] are 
for elliptic modulus k = 0.8, implying a solution of moderate amplitude. The comparison 
of these two figures serves to add credence to the relevance of the results obtained using 
the zero-amplitude solutions at the start of the bifurcation branch. 


5. With n + jjL = k and m -|- /x = /c -|- Z, for some Z G Z, the collision condition (33) is written 
as 


-F 3kl + 3k‘^ -1 = 0, 


(50) 
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(a) (b) 

Figure 6: (a) The profile of a 27r-periodic small-amplitude traveling wave solution of the 

Whitham equation Q with c « 0.7697166847, computed using a cosine collocation method 
with 128 Fourier modes, see [H]. (b) The stability spectrum of this solution, computed using 
the Fourier-Floquet-Hill method m with 128 modes and 2000 different values of the Floquet 
parameter /r. The presence of a modulational instability is clear, but no high-frequency insta¬ 
bilities are observed, in agreement with the theory presented. Note that the hallmark bubbles 
of instability were looked for far outside of the region displayed here. 


where the trivial solution I = 0 has been discarded. This is the equation of an ellipsoid 
in the {k,l) plane. It intersects lines of nonzero integer I in six integer points: ±(1,—2), 
±(0,1), ±(1, —1). Since for all of these, Q{k) = 0, any collisions happen only at the origin 
= 0. This is also illustrated in Fig. [^. 

6. The final step of our process is preempted by the results of the previous step. No Krein 
signature of colliding eigenvalues can be computed, since no eigenvalues collide. 

Since eigenvalues do not collide away from the origin they cannot leave the imaginary axis 
through such collisions and no high-frequency instabilities occur for small amplitude solutions 
of the gKdV equation. This result applies to the KdV and inKdV equations as special cases. 
The absence of high-frequency instabilities for small amplitude solutions is consistent with the 
results in, for instance, UMM- 


3.2 The Whitham equation 

As our second scalar example, we consider the Whitham equation Q. For this example, no 
analytical results exist, but the work of Sanford et al. m allows for a comparison with numerical 
results. Sanford et al. do not report the presence of high-frequency instabilities for solutions 
of any period. Their absence has been verified by us using the same methods, see Fig. Hur 
&: Johnson m consider periodic solutions, focusing on the modulational instability. However, 
they do include a Krein signature calculation of the eigenvalues of the zero-amplitude solutions, 
reaching the same conclusions obtained below. In what follows the nonlinear term N{u) does 
not contribute, as in the previous examples. 


1. The modified Hamiltonian is 




r2n 


u^dx 



y)u{x)u{y)dxdy. 
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Figure 7: (a) The dispersion relation for the Whitham equation (curve), together with the 

line through the origin of slope a;(l)/l, representing the right-hand side of ( |33[ ). (b) The curves 
-|- re), for various (integer) values of re, illustrating that collisions occur at the origin only. 


We use V to denote the speed of the traveling wave, to avoid confusion with the phase 
speed c(k) in the kernel of the Whitham equation. 

2. The dispersion relation is given by co(k) = sign(k)gk tanh(kh). 

3. The bifurcation branch starts at (17,0) = (Y^( 7 tanh(/i),0), where we have chosen N = 1 
so that the minimal period of the solutions is 27r. 

4. The elements of the stability spectrum are given by 


+ Ai)\/ 5 tanh(/i) — isign(re -|- lJ^)\/g{n + g) tanh(/i(re + g)). (51) 

5. The dispersion relation for the Whitham equation is plotted in Fig. [^a), together with 
the line through the origin with slope a;(l)/l. Since the dispersion relation is concave 


down (up) in the first (third) quadrant, the condition (33) is not satisfied. Thus collisions 
of eigenvalues away from the origin do not occur. This is also illustrated in Fig. Rb), 
where the imaginary part of An is plotted for various integers re. 

6. No Krein signature calculation is relevant since eigenvalues do not collide away from the 
origin. 


We conclude that periodic solutions of sufficiently small amplitude of the Whitham equation 
are not susceptible to high-frequency instabilities. This is consistent with the results presented 
in m, see also Fig. Thus, the Whitham equation is unable to replicate the instabilities 
found in the shallow depth water wave problem for solutions of small amplitude, despite having 
a dispersion relation that is identical to one branch of the water wave dispersion relation. We 
return to this in the next section. 


4 Two-component Hamiltonian PDEs with canonical Poisson 

structure 

We generalize the ideas of the previous section to the setting of two-component Hamiltonian 
PDEs with canonical Poisson structure. In other words, the evolution PDE can be written as 
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q 

dt \ p 


JVH ^ 



6H 



6q 


(52) 


where the Poisson operator J is given by ([^. This Poisson operator is nonsingular, thus there 
are no Casimirs. Examples of systems of this form are the Nonlinear Schrodinger equation in 
real coordinates [7], the Sine-Gordon equation [2I1ES1E0], and the water wave problem m- 
As before, our interest is in the linearization of this system around the zero-amplitude solution. 
The quadratic Hamiltonian corresponding to this linearization can be written as 



' oo ^ oo oo 

^ X] CnqL +nYl ^nPlx+p'^ 


(^riQnx 1 dx^ 


n=0 


n=0 


n=0 


(53) 


with an, bn, Cn G M. Typically the number of terms in the sums above is finite, but an example 
like the water wave problem requires the possibility of an infinite number of nonzero contributing 
terms in the Hamiltonian. As for the Whitham equation, convergence of the resulting series is 
not problematic. The form ( |53[ ) is the most general form of a quadratic Hamiltonian depending 
on two functions q{x,t) and p{x,t). Indeed, any quadratic term of a form not included above 
is reduced to a term that is included by straightforward integration by parts. The linearization 
of (52) is given by 


OO OO 

^ ^ O-nQnx T ^ ^ ( 1) bnP2nx, 

n=0 n=0 

oo oo 

Pt — ^ ^ ( 1 ) C.nq2nx ^ ^ ( 1 ) 0,nPnx' 

n=0 n=0 

We proceed with the six step program outlined in the introduction. 

1. Quadratic Hamiltonian. The modified Hamiltonian is given by 


(54a) 

(54b) 


H^ = 


^2tt / ^ oo ^ oo oo \ 

/ cpq^ 2 ^ 2 ^ ~^pY^ 

\ n=0 n=0 n=0 / 


(55) 


This expression serves as a repository for the coefficients which are needed in what follows. 


2. Dispersion relation. We look for solutions to (54a) of the form q = qexp{ikx — ioot), 
p = pexp{ikx — iujt). Requiring the existence of non-trivial {i.e., non-zero) solutions, we 
find that u}{k) is determined by 


det 


OO 


iuj + > aniiky 


n=0 


oo 

E 

n=0 


bnk 


2n 


\ 




= 0 . 


(56) 


n=0 


n=0 




This is a quadratic equation for uj{k), resulting in two branches of the dispersion relation, 
uji{k) and uj 2 {k). Assuming that ( |54a -b) is dispersive, 0 Ji{k) and 0 J 2 {k) are real-valued for 
A: G M. This is not easily translated in a condition on the coefficients an, bn, Cn and dn, 
since their reality is not assumed. 
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3. Bifurcation branches. Traveling wave solutions are stationary solutions of 


OO OO 

Qt — (^Qx T ^ ^ OjnQnx T ^ 1) bnP2nx 

n=0 n=0 


5p ’ 


OO OO 

Pt — CPx ^ 1) CnQ2nx ^ 1) 0,nPnx — 

n=0 n=0 


6J^ 

6q • 


(57a) 

(57b) 


This system has the dispersion relations Qi^ 2 {k) = uji^ 2 {k) — ck. In Fourier space the 
stationary equations become 


0 = ikcq + ^ an{ik"-)q + '^{-l)'^bn{ik)‘^"‘p, (58a) 

n=0 n=0 

OO OO 

0 = ikcp — ^(—l)"'Cn(iA:)^"'g — ^^{—l)^aniik)"^p. (58b) 

n=0 n=0 


Thus c is obtained from the condition that these equations have a nontrivial solution {q,p). 
This condition requires that the 2x2 determinant of the system above is zero. A simple 
comparison with ( |56[ ) gives that there are two bifurcation points given by {oJi{N)/N,0) 
and {uj 2 {N)/N, 0). Any positive integer value of N is allowed, but we usually choose N = 1 
so that the fundamental period is 27r. In what follows, we examine the small-amplitude 
solutions starting from the branch (c, 0) = (u;i(A')/iV, 0), without loss of generality. 


In many systems the two bifurcation branches are reflections of each other about the 
vertical axis. The corresponding solution profiles are identical to each other, moving to 
the right on one branch, moving to the left on the other. Examples are given below. Non- 
symmetric bifurcation branches cannot be excluded, however, without imposing extra 
assumptions on the coefficients of (54a-b). 


4. Stability spectrum. To find the stability spectrum, we let q = Q(x)exp(At), p = 
P(x)exp(Af). Next, using Floquet’s Theorem, 




,l]X 


(59) 


j=-oo 


j=-oo 


with p G (—1/2,1/2]. Since (57a-b) has constant coefficients, it suffices to consider 
monochromatic waves, i.e., only one term of the sums in (59) is retained. It follows 


that A satisfies (56) with ioj replaced by —A -|- i{n + p)c. Thus 


= i{n +p)c-iui{n +p) =-iQi{n +p), 1 = 1,2. (60) 

As expected, the zero solution is neutrally stable since u}ip{k) are real, assuming dispersive 
equations. The stability spectrum consists of two one-parameter point sets, one for I = 1, 
the other for 1 = 2. 

5. Collision conditions. Ignoring collisions at the origin, we require ^ 0 for 

some ni, 71-2 € 1^, p € (—1/2,1/2], li, I 2 G {1, 2}. This gives 
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Figure 8 : The graphical interpretation of the collision condition (61). The dashed curves are 
the graphs of the dispersion relations and uj 2 {k). The slope of the segment P 1 P 2 is the 

right-hand side in (61). The collision condition (61) seeks points whose abscissas are an integer 
apart, so that at least one of the segments T 3 T 4 , P^Pq-, P 5 P 4 or P^Pq is parallel to the segment 
F1P2. 


uii^jni + fi) - uji^{n2 + IJ.) ^ cJi(jV) 
ni — n2 N ' 

The right-hand side depends on cji since we have chosen the hrst branch of the dispersion 
relation in Step 3. As before, this collision condition may be interpreted as a parallel 
secant condition, but with the additional freedom of being able to use points from both 
branches of the dispersion relation. This is illustrated in Fig. 

6. Krein signature. In the setting of a system of Hamiltonian PDFs as opposed to a scalar 
PDF, we use a different but equivalent characterization of the Krein signature |35]. The 
Krein signature is the contribution to the Hamiltonian of the mode involved with the 
collision. Since our Hamiltonians are quadratic, this implies that the Krein signature of 
the eigenvalue A with eigenvector v is given by 


signature(A, u) = sign(u^£c^^)) 


(62) 


where Cc is the Hessian of the Hamiltonian P[^, and denotes the complex conjugate of 
the transposed vector. Since the Hessian Cc is a symmetric linear operator, the argument 
of the sign in (62) is real. Recall that the linearization of the system (§ can be written as 


dt 


= JCc 


( 63 ) 


which makes it easy to read off Cc- For the case of (57a-b), 
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(64) 


/ 


/:c = 


V 


Y.Cn{-lTdl 

n=0 

oo 

cdx + ^ and. 


-cdx + '^an{-lTdl 

n=0 

oo 

Y^hni-lTdl 


^2n 

■^x 


n=0 


n=0 


Next, the eigenvectors v are given by 


^ \ _ \t-\-i^x-\-inx I '^n 

p j \ P. 


Qr 


where {Qn,Pn)'^ satisfies 


f I = JCc 


Qn 

Pn 


(64): 


(65) 


( 66 ) 


Here Cc is the symbol of £c, *-e., the 2x2 matrix obtained by replacing dx —)■ i{n + p) in 


Cr = 


'^Cn{n + 


-ic{n + p) + '^ an{-l)^iin + ipY 


n=0 


n=0 


ic{n + p) + '^ an{in + ipY 


'^bnin + pY"^ 


■ (67) 


n=0 


n=0 


Using (66), (62) is rewritten as 


signature = sign det 


Qn Pn 

Q * p* 

n 


( 68 ) 


The determinant is imaginary, since interchanging the rows gives the complex conjugate 
result. Since i® imaginary, the result is real and the signature is well defined. Again, 

it is clear that no conclusions can be drawn if Y'Y = 0. Since we wish to examine whether 

n,l 

signatures are equal or opposite, we consider the product of the signatures corresponding 
to ■ Using (67) we find that signatures are opposite, provided that 

/t-l 7t2 5^2 


OO 


OO 


N + pf^^ Y1 + P?^"" 

ii=o i2=o 


^lYni + p) + 02 ^ 3+1 (-l)^Hni +/r)2^3+ij 


X 


^p{n2 + p) + Y 0'-2j4+ii-'^yYn2 + 


< 0 . 


(69) 


The above condition is obtained by expressing the eigenvectors in (66) in terms of the 
entries of the first row of (67). An equivalent condition is obtained using the second row: 
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oo 


oo 


bj.im + Y 

ji=0 j2=0 


j Ui^im + fi) - Y «2i3+i(-lF®(^i ] 

V i3=o / 


X 


^l2{n2 + tJ-)-Y 02 ^ 4+1 (- 1)-^^(^2 + < 0 . 


(70) 


Depending on the system at hand, the condition (69) or (70) may be more convenient to 
use. 


Remark. An important class of systems is those for which uji{k) = —uj 2 {k). We refer to 


such systems as even systems. It follows immediately from (56) that for even systems 02^+1 = 0, 


j = 1,2,.... The Krein conditions (69) and (70) simplify significantly, becoming 


w«i(ni + lj)u:i^{n2 + k-) Y Y '^^ 2(^2 + kf^^ < 0, 

jl=0 j2=0 


or 


(71) 


{ni + k)^i2{^2 + k)Y < 0- (72) 

jl=0 j2=0 

We summarize our results. 


Assume that the linearization of the Hamiltonian system (52) is dispersive {i.e., its disper¬ 
sion relations u)i{k) and 0 J 2 {k) are real valued for A: G M). Let iV be a strictly positive integer. 
Consider 27r/A^-periodic traveling wave solutions of this system of sufficiently small-amplitude 
and with velocity sufficiently close to a;i(A^) /N. In order for these solutions to be spectrally un¬ 
stable with respect to high-frequency instabilities as a consequence of two-eigenvalue collisions, 
it is necessary that there exist h, I 2 G {1, 2}, ni, n 2 G Z, ni 7 ^ re 2 , /r G (—1/2,1/2] for which 


wzi (ni + n) u!i (N) u}i2 (n2 + k) 

/ AT ’ I / 


ni + 


N 


n2 + k 


N 


such that 


(^1 + k) - (^2 + k) _ ^i(iv) 

ni — n2 N 


(73) 


(74) 


and (69), or equivalently, (70) holds. 


We proceed with examples. 


4.1 The Sine-Gordon equation 

As a first example, we consider the Sine-Gordon (SG) equation |32] : 


utt - Uxx + smu = 0. (75) 

The stability of the periodic traveling wave solutions of this equation has been studied recently 
by Jones et al. [23 EO]. Different classes of periodic traveling wave solutions exist, but only two 
of those can be considered as small-amplitude perturbations of a constant background state. 
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We consider the so-called superluminal (c^ > 1) librational waves. The subluminal (c^ < 1) 
librational waves require the use of the transformation v = u — n so that their small amplitude 
limit approaches the zero solution. We do not consider them here. The limits of the rotational 
waves are either soliton solutions or have increasingly larger amplitude. As such the rotational 
waves do not fit in the framework of this paper. An overview of the properties of these solutions 
as well as illuminating phase-plane plots are found in [2^. In contrast to 123 EO], we fix the 
period of our solutions, as elsewhere in this paper. This makes a comparison of the results more 
complicated. 


1. Quadratic Hamiltonian. With q = u, p = ut. 


r2n 


H^ = 


cpQx + 1 dx. 


Thus 6o = 1) Co = 1, Cl = 1 are the only non-zero coefficients. 


2. Dispersion relation. Using (56), 


(76) 


Wi,2 = ±Vl + A:2. (77) 

These expressions are real valued for k £ M., thus the SG equation is dispersive when 
linearized around the superluminal librational waves. Both branches of the dispersion 
relation are displayed in Fig. &• 

3. Bifurcation branches. With = 1, we obtain c = wi(l)/l = ^/ 2 . 

4. Stability spectrum. The stability spectrum is given by (|60[): 


= -iVli{n + p) = i{n + p)V2 =F i^/l + {n + pY, (78) 

with / = 1 (/ = 2) corresponding to the — (-|-) sign. Here n £ 7j, p £ [—1/2, 1/2). 

5. Collision condition. The collision condition ( |61[ ) becomes 

ui{ni + p) - u]2{n2 + p) ^ ^ 

ni - 712 

We have chosen wq = oji and = a; 2 , since it is clear that the collision condition can only 
be satisfied if points from both dispersion relation branches are used. This is illustrated 
in Fig.|^. In fact, many collisions occur, as is illustrated in Fig. [33. One explicit solution 
is given by 


Til = 3, ?12 = 0, p = 


VlO- 3 


0.081138830. 


(80) 


6. Krein signature. Since uj 2 {k) = —u}i{k), we may use the conditions © or ( [T^ )- Since 
only one bj ^ 0, (72) is (slightly) simpler to use. We get that 


iOi{ni + p)uj2{n2 + p) < 0 (81) 

is a necessary condition for the presence of high-frequency instabilities of small-amplitude 
superluminal librational solutions of the SG equation. The condition is trivially satisfied 
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(b) 


Figure 9: (a) The two branches of the dispersion relation for the Sine-Gordon equation. The 

line segment P 1 P 2 has slope t<;(l)/l, representing the right-hand side of (61). The slope of the 
parallel line segment P 3 P 4 represents the left-hand side of (61). (b) The two families of curves 
Q,i{k+n) (red, solid) and Q 2 {k + n) (black, dashed), for various (integer) values of n, illustrating 
that many collisions occur away from the origin. 


as it was remarked in the previous step that points from both dispersion relation branches 
have to be used to have collisions. 


It follows that for the superluminal solutions of the SG equations the necessary condition for 
the occurrence of high-frequency instabilities is satisfied. Nevertheless, as the results of 


show, such instabilities do not occur. This is illustrated in Fig. 10 The left panel illustrates an 
exact 27r-periodic superluminal solution of the SG equation with c ss 1.236084655663, obtained 
using elliptic functions. Gomputing the stability spectrum (right panel) of the solution using the 
Fourier-Floquet-Hill method m with 51 Fourier modes and 1000 Floquet exponents shows that 
no high-frequency instabilities are present to within the accuracy of the numerical method. This 
is consistent with the results of [291130] where only the presence of a modulational instability is 
observed. Thus the example of this section illustrates that the necessary condition is not always 
sufficient. 


4.2 The water wave problem 

As a final example, we consider the water wave problem: the problem of determining the 
dynamics of the surface of an incompressible, irrotational fluid under the influence of gravity. 
For this example, the effects of surface tension are ignored and we consider only two-dimensional 
fluids, i.e., the surface is one dimensional. The Euler equations governing the dynamics are 


4^xx 4^zz — 0? 

(x,z) G D, 

(82a) 

II 

z = —h, 

(82b) 

Vt Vx4^X — ^z^ 

II 

(82c) 

+ <Pz) + gv = 0 , 

II 

(82d) 

where x and z are the horizontal and vertical coordinate, respectively, see Fig. 11 

z = r]{x, t) is 


the free top boundary and 4){x, z, t) is the velocity potential. Further, g is the acceleration due 
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(a) (b) 

Figure 10: (a) A small-amplitude 27r-periodic superluminal solution of the SG equation 

(c 1.236084655663). (b) A blow-up of the numerically computed stability spectrum in a 
neighborhood of the origin, illustrating the presence of a modulational instability, but the ab¬ 
sence of high-frequency instabilities. 





Figure 11; The domain for the water wave problem. Here z = 0 is the equation of the surface 
for flat water, z = —h is the flat bottom. 


to gravity and h is the average depth of the fluid. 

The main goal of the water wave problem is to understand the dynamics of the free surface 
r](x,t). Thus it is convenient to recast the problem so as to involve only surface variables. 
Zakharov m showed that the water wave problem is Hamiltonian with canonical variables 
r]{x,t) and (p{x,t) = 4>{x,r]{x,t),t). In other words ip{x,t) is the velocity potential evaluated at 
the surface. Following [12], the Hamiltonian is written as 

1 

^ = 2 Jq 

where G{r]) is the Dirichlet —)• Neumann operator: G{r])(p = (1 -|-at z = r]{x,t). Here 
(pn is the normal derivative of p. Using the water wave problem, G{r])ip = p^ — rjxPx = Vt: which 
is the first of Hamilton’s equations. The water wave problem for ri{x, t) and (p{x, t) is 


Vt = 


5H_ 

6(p ’ 


= 


6 r] ' 


(84) 
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1. Quadratic Hamiltonian. Since for our purposes, the linearization of (84 b) in a moving 
frame is required, it suffices to evaluate the Dirichlet —)• Neumann operator G{r]) at the 
flat surface i] = 0, resulting in 


G{0) = —idx tanh.{—ihdx)■ 
The quadratic Hamiltonian is given by 


(85) 


hO = c 


I'27T I'27V 

I ipr]xdx + - {(p{-itanh.{-ihdx)ipx + gv'^) dx, 
0 ^ Jo 


( 86 ) 


giving rise to the linearized equations in a frame moving with velocity c: 


rjt = —— = crjx - itanh[-ihdx)(px, 


= 


6cp 

51^ 

6r] 


= apx + gr]. 


(87a) 

(87b) 


2. Dispersion relation. The well-known dispersion relation [4^ for the water wave problem 
is immediately recovered from the linearized system (87a-b) with c = 0 (no moving frame), 
resulting in 


ui‘^ = gktanh.{kh). (88) 

Note that the right-hand side of this expression is always positive. Thus there are two 
branches to the dispersion relation: 


coi ^2 = ±sign(A:) -\/gk tanh(A:/i). (89) 

Thus oji {L 02 ) corresponds to positive (negative) phase speed, independent of the sign of 
k. 


3. 


4. 


Bifurcation branches. Branches originate from (c, amplitude) = (a;i(l)/l,0) and 
(c, amplitude) = (a; 2 (l)/l, 0). Without loss of generality, we focus on the hrst branch 
for which the phase speed gtanh{h) is positive. This allows for a straightforward com¬ 


parison of our results with those for the Whitham equation, in Example 3.2 
Stability spectrum. The elements of the spectrum are given by 


= -i^i{n + g) 

= i(n-|-/i)y^5Tanh(/)^ — isign(n-|-(90a) 
= -i^2{n + n) 

= i{n + g) \/g tan]i{h) + isign(n -|- fj,)\/g{n + g) tanh(/i(n -|- g)). (90b) 


The sign(n-|-^)’s may be omitted in these expressions, as the same set of spectral elements 
is obtained. 
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(c) 


Figure 12: (a) The two branches of the dispersion relation for the water wave problem [g = 1, 

h = l). The line through the origin has slope t<;i(l)/l, representing the right-hand side of (61). 
(b) The two families of curves Hi (A: -|- n) (red, solid) and Vt. 2 {k + n) (black, dashed), for various 
(integer) values of n, illustrating that many collisions occur away from the origin, (c) The origin 
of the high-frequency instability closest to the origin as a function of depth h. 


6 . 


5. Collision condition. The condition (61) is easily written out explicitly, but for our 


purposes it suffices to plot Qi{k + n) and il. 2 {k + n), for different values of n. This is done 
in Fig. & with g = 1 and h = 1. Although only the first collision is visible in the figure 
(all intersection points are horizontal integer shift of each other and correspond to the 
same value of g and it is clear from the curves shown that many collisions occur. 

The figure is qualitatively the same for all finite values of depth. 

Krein signature. The conditions 0 and ([72|) become 


ug{ni + g)uji^{n2 + g)g^ < 0 , 


(91) 


and 


wq(ni -F g)uji^{n2 + Z^) ^ ^ ^(n2 -h gf^^ < 0, (92) 

jl = l j2 = l 
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respectively. Here the coefficients aj are related to the Bernoulli numbers m. as they 
are defined by the Taylor series 


tanh( 2 ;) = , \z\ < 7r/2. (93) 

j=0 


Because of the finite radius of convergence of this series, (92) is only valid for small values 
of the wave numbers ni + /r and n 2 + but it is possible to phrase all results in terms of 
tanh directly, avoiding this difficulty. For instance, using (93), (92) may be rewritten as 


+ ^)uji^{n2 + //) 


Wa(ni + fi) ^‘lin2 + 


< 0 

9 9 

i^h (^1 + (”2 + /r) < 0, 


in agreement with ( |91[ ). The indices a, /3 G {0,1} are irrelevant since uja and both 
appear squared. This serves to illustrate that for specific examples one of the two criteria 
(69) and (70) (or ( [7l| ) and (72) for even systems) may be significantly easier to evaluate, 
although they are equivalent. 


Thus all collision points are potential origins of high-frequency instabilities. It appears 
from the numerical results in m that the bubble of non-imaginary eigenvalues closest 
to the origin contains the high-frequency eigenvalues with the largest real part. Thus 
for waves in shallow water kh < 1.363 (no Benjamin-Feir instability) [U WI\ 154] . these 
are the dominant instabilities. For waves in deep water (kh > 1.363) the Benjamin-Feir 
instability typically dominates, although there is a range of depth in deep water where the 
high-frequency instabilities have a larger growth rate, see m- The dependence on depth 
h of the location on the imaginary axis from which the high-frequency bubble closest to 
the origin bifurcates is shown in Fig. [T2)(c), with 5 = 1. As /i —?• oo, the imaginary part of 
A —)• 3/4. This asymptote is drawn in Fig. |12[c) for reference. This figure demonstrates 
that for all positive values of the depth h, the instabilities considered are not modulational 
as they do not bifurcate away from the origin as the amplitude increases. 


It was remarked in Example 3.2 that no collisions are possible due to the concavity of 
the dispersion relation. As a consequence, all collisions away from the origin observed in 
Fig.[T^ involve both branches of the dispersion relation, i.e., they involve a solid curve 
and a dashed curve. This is easily seen from Fig.|12^: a parallel cord with abscissae of the 
endpoints that are integers apart is easily found by sliding a parallel cord away from the 
cord ((0,0), (1, a;i(l))) until the integer condition is met. This implies that {ni+^) and 
w/ 2 (?T '2 + k) in the collision condition (61) have opposite sign and (91) is always satisfied. 
Thus colliding eigenvalues of zero-amplitude water wave solutions always have opposite 
Krein signature. As a consequence, the necessary condition for the presence of high- 
frequency instabilities is met. In fact, it was observed in m that all colliding eigenvalues 
give rise to bubbles of instabilities as the amplitude is increased. 


Our general framework easily recovers the results of MacKay & Saffman |36j . There the 
set-up is for arbitrary amplitudes of the traveling wave solutions, but the results are only 
truly practical for the zero-amplitude case. 


Remark. It follows from these considerations that the high-frequency instabilities present in 
the water wave problem are a consequence of counter-propagating waves as no such instabilities 
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are present in the Whitham eqnation Q. Althongh it is often stated that the value of the 
Whitham equation lies in that it has the same dispersion relation as the water wave problem 
(see for instance [50]), this is in fact not the case as it contains only one branch of the dispersion 
relation. Thus the equation does not allow for the interaction of counter-propagating modes, 
and as such misses out on much of the important dynamics of the Euler equations. 


5 A Boussinesq-Whitham equation 


The goal of this section is the introduction of a model equation that has the same dispersion 


relation as the Euler equations (82a-d) at the level of heuristics that led Whitham to the 
model equation Q. In other words, we propose a bidirectional Whitham equation, so as to 
capture both branches of the water wave dispersion relation. We refer to this equation as the 
Boussinesq-Whitham (BW) equation. It is given by 


qu = N{q) + &: 


(aq^+ r K{ 

\ J —OO 


X - y)q{y)dy 


(94) 


where 


= 2; 


c\k)e^’^^dk, 


(95) 


and c^(A:) = gtanh{kh)/k, a > 0 for the water wave problem without surface tension. In (94), 


N{u) denotes the nonlinear terms, which are ignored in the remainder of this section. Since 
our methods focus on the analysis of zero amplitude solutions, the sign of a is not relevant in 
what follows. This equation is one of many that may stake its claim to the name “Boussinesq- 
Whitham equation”. Equation (94) is a “Whithamized” version of the standard Bad Boussinesq 


equation and it may be anticipated that it captures at least the small-amplitude instabilities 
of the water wave problem in shallow water. It should be remarked that the Bad Boussinesq 
equation is ill posed as an initial-value problem |3Z|, but it might be anticipated that the 
inclusion of the entire water-wave dispersion relation overcomes the unbounded growth that is 
present due to the polynomial truncation. We return to this at the end of this section. 

Before applying our method to examine the potential presence of high-frequency instabilities 
of small-amplitude solutions of the BW equation, we need to present its Hamiltonian structure. 


Eurther, since (94) is defined as an equation on the whole line, a periodic analogue is required, 
as in Section O 


It is easily verified that (94) is Hamiltonian with (non-canonical) Poisson operator [37] 


J = 


0 9 ^, 

dx 0 


(96) 


and Hamiltonian 


H = 


P^ + ^q^]dx + -J dx J dyK{x-y)q{x)q{y). (97) 


Indeed, (94) can be rewritten in the form (llj) with u = {q,pY ■ 


To dehne a periodic version of (94), let 


^ OO 

K{x) = 


(98) 


j=-oo 
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where kj = 271^1 L, j G Z. The periodic BW equation is obtained from Q, using ( [%| ) and 
(97), but with all ± infinities in the integration bounds replaced by ±L/2, respectively. Since 
(94) has a Poisson operator (96) that is different from those used in Sections and minor 
modifications to the use of the method are necessary. 

1. Quadratic Hamiltonian. Ignoring the contributions of the nonlinear term, the 
quadratic Hamiltonian in a frame of reference moving with speed V is given by 


Hy - 


f27r 


1 . 

Vqp+ -p‘ 


dx+\ 


p27T r27r 

/ dx dyK{x-y)q{x)q{y), 
'o Jo 


(99) 


where we have fixed the period of the solutions to be L = 27r. The inclusion of the first 


term in (99) is one place where the effect of the different form for J is felt, as its functional 


form is a direct consequence of the form of (96). 


2. Dispersion Relation. A direct calculation confirms that 


uj'^ = gktanh.{kh), (100) 

which is, by construction, identical to the dispersion relation for the full water wave 


problem (88). This gives rise to two branches of the dispersion relation (89), corresponding 


to right- and left-going waves. 

3. Bifurcation Branches. Bifurcation branches for 27r-periodic solutions start at (Hi,2)0), 
where the phase speeds Vi ,2 are given by Vi ,2 = ±y/ 5 tanh(/i). 

4. Stability Spectrum. The stability spectrum elements are, again by construction, iden¬ 


tical to those for the water wave problem, given in (90a-b). 


5. Collision Condition. Given that the spectral elements are identical to those for the 
water wave problem, the collision condition is identical too. It is displayed in Fig. [T^a-b). 
Thus, collisions away from the origin occur. It remains to be seen whether these can result 
in the birth of high-frequency instabilities. 

6. Krein Signature. As for the canonical case of Section]^ we use (62). Thus we calculate 
the Hessian Cr of the Hamiltonian 


Let 


j=0 


( 101 ) 


where yj = with the coefficients aj defined in (93). A direct calculation gives 


that the Hamiltonian (99) is rewritten as 


Hy - 


r2n 


p^ + Vqp + 


OO 

I dx- 

j=0 


( 102 ) 


Using this form of the Hamiltonian, the calculation of the Hessian is straightforward, 
leading to 
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(103) 


Ly = 


i=o 

V 1 


Next, we compute the eigenvectors v = {q,p)'^- We have 


9 A jxt ( Q{x) 

p ) I P{x) 


(104) 


where {Q, P)^ satisfies 


A 


Q 

P 


0 dx \ ^ [ Q 


dx 0 


Cy 


P 


(105) 


This is a second place where the Poisson operator J plays a crucial role as it affects the 
form of u = {q,p)'^ and thus the expression for the signature. One easily verihes that 


Q \ ^ Jin+i,)x ( i{n + p) 

P J \\-i{n + p)V 


satisfies (105). 


We need to evaluate the sign of 


(106) 


Q Vr ( Q 

p ) p 




-i{n + p) 

—X + i{n + p,)V 


-i{n + p) 

—X + i{n + p)V 


t( 

V 


t( 

V 


j=0 

V 1 


i{n + p) 

X — i{n + p)V 


o^n+px 


+ y \ ( i{n + p) 

i=o I I A — i{n + p)V 


V 


1 


-i(n + p) A / c^(n + p) y \ f i{n + p) 
iuj{n + p) J y V 1 J \ —iu}{n + p) 

= 2w (w — (n + p)V). 


(107) 


Let the signature associated with the first eigenvalue be the sign of 2a;j^(t(jj^ — (rei + p)V), 
where is a function of ni + p. Similarly, for the second eigenvalue, the signature 
is the sign of 2ujj^{iOj2 — (n 2 + p)y)- Using the collision condition the 

product of these two expressions is 4:Ujj^ijjj2{ujj2 — in 2 + p)yy, which is less than zero since 
collisions can only occur for eigenvalues associated with opposite branches of the dispersion 
relation, see Fig. It follows that, as in the water wave case, the signatures of colliding 
eigenvalues are always opposite, and the necessary condition for spectral instability is 
met. Thus, unlike the Whitham equation (1^, the BW model (94) does not exclude the 


presence of high-frequency instabilities of small-amplitude solutions. 
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The results obtained from the Krein signature calculations are confirmed by numerical re- 
Panel (a) shows a numerically computed traveling wave solution of the BW 


suits, see Fig. 13 


equation (94). This solution is computed using a cosine collocation method with 60 points, as 
for the Whitham equation, see Fig. ^ [H]. For the solution plotted, c ~ 1.049815. The second 


panel in the first row displays the spectrum computed using Hill’s method with 100 modes and 
20000 values of the Floquet parameter, using an interpolation of the solution profile. This panel 
shows the presence of a large number of apparent instabilities, most with small growth rate, in 
the neighborhood of the imaginary axis. The third panel shows a zoom of the region around 
the origin, revealing a modulational instability. This is expected, since such an instability is 
also present for the Whitham equation, see Section 3.2 The fourth panel zooms in on the first 


bubble of instabilities centered on the positive imaginary axis, revealing a shape and location 
that is consistent with the Krein collision theory presented here. 



X 


(a) 


ImA 



ReA 


(b) 



(c) (d) 


Figure 13: (a) A small-amplitude traveling wave solution of the Boussines-Whitham equation 

(94) with c ~ 1.0498515. (b) The numerically computer stability spectrum, (c) A blow-up of 
the stability spectrum in a neighborhood of the origin, (d) A blow-up of the stability spectrum 
around what appears as a horizontal segment visible in (b) immediately above the longest 
segment appearing horizontal. More detail is given in the main text. 


The wave form displayed in Fig. [l^a) does not have zero average, unlike the one shown 
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in Fig. for reasons explained here. Let us examine the stability of a flat-water state q = a 
(constant), p = 0. Thinking of the BW equation as an approximation to the water wave problem 
where the flat-water state is neutrally stable (spectrum on the imaginary axis), independent of 
the reference level of the water, the neutral stability of this state is desired in the context of the 
BW system as well. However, the BW system is easily checked to not be Galilean invariant, 
thus the average value of the solution may be important. 

Linearizing the system around the flat water state {q, p) = (a, 0) results in a linear system 
with constant coefficients, whose dispersion relation is given by = 2ak‘^+k'^c^{k). This results 
in two branches for the dispersion relation: uii ^2 = ^k^yc^{k) + 2a. It follows that if o > 0 then 
both oji and ijJ 2 are real, resulting in neutral stability, since the stability eigenvalue and the 
frequency uj{k) are different by a factor of i. On the other hand, if a < 0, it follows that both oji 
and 0 J 2 are imaginary for sufficiently large k, since lim|fc|_,.oo c{k) = 0. This leads to the dynamics 
of the flat-water state with a < 0 to not only be unstable, but to be ill-posed, as the growth 
rate of the instability —>■ oo as |fc| —)> oo. Thus Whitham-izing the Bad Boussinesq equation and 
incorporating the full water wave dispersion relation does not remove the illposedness of the 
problem. Rather it alters it where negative constant solutions experience unbounded growth, 
unlike positive constant solutions. 

It is observed numerically that this behavior of perturbed constant solutions is carried 
over to nonconstant solutions: solutions of negative average display the same illposed behavior 
described above, with stability spectra that have unbounded real part. In contrast, the spectra 
of solutions of positive average have bounded real part, as in Fig. 13 Annoyingly, the illposed 


behavior extends to numerical solutions constructed to have zero average. Presumably this is 
a consequence of numerical error, as higher accuracy numerical experiments display narrower 
spectra whose real part tends to infinity more slowly. 

We may summarize our findings on the BW equation as follows. The equation was con¬ 
structed as a bi-directional Whitham equation so as to truly have the same linear dispersion 
relation as the water wave problem. Even though the BW has a different Poisson structure 
than the water wave problem, we find that periodic solutions of the BW are susceptible to 
high-frequency instabilities originating from Krein collisions at the exact same locations on the 
imaginary axis as the water wave problem. On the other hand, we have not attempted to quan¬ 
tify whether the resulting growth rates are comparable to those for the water wave problem. 
Further, the illposedness of the equation for solutions of negative average is a significant strike 
against its potential use in applications. Nevertheless, it appears possible to design more equa¬ 
tions like the BW equation, possessing the exact same dispersion relation as the water wave 
problem and with it all its high-frequency instabilities, without the equation dynamics being 
illposed for any important class of solutions. 
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